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INTRODUCTION 

The accurate calculation of thermal stresses in complicated airframe 
structures often requires a refined finite-element grid and a corresponding 
large computational time. In looking toward combined thermal -structural 
design and optimization calculations, it is essential to avoid having to solve 
an excessively large system of equations since such calculations are performed 
many times during optimization procedures. 

A procedure for generating accurate thermal stresses with coarse 
finite-element grids was developed and described by Ojalvo (ref. 1). The 
procedure is based on the observation that for linear thermoelastic problems, 
the thermal stresses may be envisioned as being composed of two 
contributions--the first due to the strains in the structure which depend on 
the integral of the temperature distribution over the finite element and the 
second due to the local variation of the temperature in the element. Ojalvo's 
key idea was that the first contribution could be accurately predicted with a 
coarse finite-element mesh. The resulting strain distribution could then be 
combined via the constitutive relations with detailed temperatures from a 
separate thermal analysis. The result would be accurate thermal stresses from 
coarse finite-element structural models even where the temperature 
distributions have sharp variations. 

Although this intriguing idea was proposed in 1974, its use has not been 
documented in the open literature except for the original AIAA Technical 
Note, It has recently received attention by this author because of the 
current interest at Langley in rapid analysis and design-oriented analysis 
techniques. The range of applicability of the method for various classes of 
thermostructural problems such as in-plane or bending type problems and the 
effect of the nature of the temperature distribution and edge constraints was 
not fully documented. These questions are addressed in this paper. Ojalvo's 
method is used in conjunction with the SPAR finite-element program (ref. 2) 
and extensive calculations are carried out and are described. Results are 
obtained for rods, membranes, a box beam and a stiffened panel. 
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SUMMARY OF OJALVO' S METHOD 


The following is a brief summary of "Ojalvo's" method (ref. 1) for thermal 
stress calculations using a finite-element (F.E.) analysis procedure. For a 
detailed explanation of the following subject topics see references 1 and 3. 
The equations shown in figure 1 represent the system of matrix equations 
involved in a typical structural F.E. analysis. The system stiffness matrix 
[K] relates the vector of element node point displacements {6} and mechanical 
and thermal load vectors ({ Pmech} {Pf} respectively) by the first 
equation. The stiffness matrix can be calculated from the strain shape 
function matrix [B] and the constitutive relations [D] as shown. If a 
structure, idealized by finite elements, has sufficient nodal constraints so 
the stiffness matrix [K] is nonsingular, the node-point displacements {6} can 
be obtained from the first equation. These displacements are then used to 
calculate corresponding strains {e} and stresses {a}. 

As explained in reference 1, the dependence of {e| upon the temperature T for 
a linear elastic analysis is through {P-|-} and {P-j-} is a function of the 
spatial integral of temperature throughout the element. Hence it appears {e} 
is related to integrals of T. The stresses {o}, however, are not only 
related to integrals of T through {e} but are also directly related to the 
local temperature. Thus, one can expect a greater accuracy in the numerical 
calculation of strain than stress in thermomechanical problems since some of 
the errors of approximations in T are self-cancelling in the integrals which 
determine strain. In essence one should calculate accurate strains with as 
coarse a F.E. grid as possible to approximate the thermal load vector. 
Theoretically one could then use a coarser structural F.E. grid than a 
corresponding thermal F.E. grid to obtain the same desired accuracies in each 
analysis. One could improve the coarse-grid F.E. stress results by using the 
coarse-grid structural results for |e} and coupling them with the fine-grid 
thermal results for T. This last statement and the last equation of the 
figure summarize Ojalvo's method. By comparison, in the conventional F.E. 
procedure there is no separation in the strain and stress calculation; the 
same F.E. grid is used to calculate consistent nodal thermal load, and nodal 
strains and stresses. Also, in most finite elements the temperature field in 
the element is similar to the displacement field; hence, sharp variations in 
temperature are not accounted for. 

FOR FINITE ELEMENT ANALYSIS: 

(6} = (Pmech) ^ (’’t) ' ,^T(x,y,z)[B:^ CD]{a} dV 

M = /Cb]^ Cd^Cb] dV 

(0 ' + ^T(x,y,z)CB]^CD]{a} dv)j 

(o) = [D] ({e} - T{a}) 

• (e) IS RELATED TO INTEGRAL OF T 

• (a) IS REUTED TO THE LOCAL TEMPERATURE 

• THEREFORE {e} SHOULD BE NUMERICALLY MORE ACCURATE THAN { 0 } 

HENCE CALCULATE STRESS AS: (a) = [O] (fE) - T/a}) 

COARSE-Cm ^FINE-GRID 
STRUCT. ANALYSIS THERMAL ANALYSIS 

Figure 1 


116 



SHUTTLE ORBITER VERTICAL FIN STRUCTURAL IDEALIZATION 
ANALYZED BY OJALVO 


To demonstrate the usefulness of Ojalvo's method, reference 1 presents results 
for a hot structure concept for the vertical fin of the Space Shuttle Orbiter 
(ref. 4). Temperatures were calculated by a detailed lumped parameter 
analysis and used as input to a coarse-grid F.E. structural analysis. As 
shown in figure 2, only a section of the fin main structural box between ribs 
9 and 10 was idealized by finite elements using rod, membrane and shear 
elements. Rene' 41 was as the structural material at locations where the 
temperature exceeded 923K (1200°F) and Inconel 718 was used for locations 
below that temperature. The spars were corrugated to partially alleviate 
excessive thermal stresses. 




STRUCTURAL IDEALIZATION OF FIN 
BETWEEN RIBS 9 AND 10 


Figure 2 
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RESULTS OF OJALVO'S METHOD FOR THE SHUTTLE ORBITER VERTICAL FIN 


A chordwise temperature distribution from a detailed lumped-parameter 
transient thermal analysis is shown as a solid curve in the upper left-hand 
part of figure 3. The massive spars and spar caps act as heat sinks causing 
sharp temperature drops in their vicinity during transient heating. The 
dashed lines are the integrated average values of temperature which were used 
as input for the cover panels in the structural F.E. model. Since the fin is 
long and slender and since the section analyzed is far enough away from the 
ends beam theory may be used. If z is the coordinate direction through the 
depth of the beam (along the chord in this example), and assuming T = T(z) and 
beam symmetry in two directions, the following equation can be used to 
predict strains (ref. 1): 


e 


/ EaTdA 
= A + 
jEdA 
A 


z/ EaTzdA 

/Ez^dA 

A 


Because of the above relationship, the sharp drops in temperature near the 
spars have little effect on the strains. Hence, as shown in figure 3, the 
strain distribution in the covers predicted by coarse-grid F.E. is linear and 
very close to the predicted beam theory solution regardless of the nature of 
the temperature distribution. Thus a coarse grid is sufficient to obtain 
accurate strains for a thermal stress problem which would require a fine grid 
to predict the temperature distribution. Also, mesh refinement of the 
structural grid may not be necessary for performing quasi -static 
thermal /structural analysis. 

A comparison of stress distributions using beam theory, conventional F.E., and 
Ojalvo's method is shown in figure 3 at the right. The conventional F.E. 
procedure, using constant-strain constant-temperature quadrilateral membrane 
elements and integrated average temperature values miss large peak stresses in 
the covers. However, stress results using the strains from the coarse-grid 
structural F.E. analysis and the existing detailed temperature distribution 
(Ojalvo's method) agree closely with beam theory results. Ojalvo's method 
works well for this problem because the strain distribution is linear and the 
thermal load vector can be closely approximated by a coarse grid (4 elements). 



Figure 3 
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QUALIFICATION OF OJALVO'S METHOD 


With regard to earlier statements in reference 1, it is true that the node 
point displacements {6} are directly related to the spatial integral of the 
temperature distribution as shown by the first equation in figure 4. However, 
the strains are directly related to the strain shape function matrix [B] which 
relates strains to node point displacements {6}, As shown in the figure [B] 
is related to the strain operator matrix [L] which performs the operation of 
differentiation and the shape function matrix [N] which relates the 
displacement vector {u} to the vector of node point displacements {6}. The 
[B] matrix operates on the {6} and transforms it into {e}; embedded in this 
operation is the differentiation of the shape function matrix [N]. The 
direction of differentiation (implied in [L]) is dependent upon the type of 
element. Hence, depending on the arguments of integration in {Pj} and the 
direction of differentiation in [L] the strains may be directly related to the 
local temperatures and not their integrals as stated in reference 1. 


As stated in reference 5, deflections caused by heating were "acceptable" but 
"special attention" was necessary for the "interpolation of the stresses which 
could deviate considerably from the true stress state." Here, reference 5 
indicates that deflections were less sensitive to temperature distributions 
than their corresponding stresses, as stated in reference 1; however, no 
mention is made of the sensitivity of strains to temperature distributions. 

It appears that for bending problems, beams, or plates, there is a direct 
relationship between e and integrals of T. The reason for this 
relationship is that for beam-type problems the argument of integration for 
|Pj} is through the depth of the beam while the desired stresses are in the 
longitudinal direction. This means the direction of differentiation implied 
by [B] is along the length and hence does not affect the integration. The 
direction of differentiation for membrane-type problems, however, is the same 
as the arguments of integration in { Pj} and hence the strain variation 
should be similar to the temperature variation. This reasoning suggests that 
while there appears to be an advantage in using Ojalvo's method for bending 
problems, it does not appear to be appropriate for membrane or plane-stress 
problems. 

(61 ■ +/T(x,y,z([B]^[0]|aldv) 

(£| ■ [B]|6I ■ :L][N](6| 

WHERE FOR EXAMPLE 


CL]{u} = 



d/dx 0 
0 d/by 
b/by d/dx. 




[L] = STRAIN OPERATOR (DIFFERENTIATION) 
MATRIX 


CN]{6| “ {u} ; N = SHAPE FUNCTION MATRIX 
HENCE [B] PERFORMS DIFFERENTIATION 

• BENDING PROBLEM Jd ~ /Tdv 

• MEMBRANE OR PLANE STRESS PROBLEM {£( ~ T 


Figure 4 
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COMPARISON OF STRESS CALCULATION METHODS FOR CLAMPED ROD PROBLEM 


To qualify the usefulness of Ojalvo's method, analytical and numerical results 
for several classes of thermostructural problems are presented beginning with 
the simplest one-dimensional example, a rod. As shown in figure 5, an 
aluminum rod of length i is clamped at both ends and subjected to a 
sinusoidal temperature variation along its length; the governing equations are 
as follows: 

( 1 ) _ a -II =0 u(o) = u{i) = 0 

d'x2 dx 

Integrating and substituting T = Tq sin ^ 

(2) u(x) = - “^0^ cos — - 

7T £ 7T 7T 

and 


(3) £x(x) = aToSin = aT - 

£. TV TV 


a = -2EaTo/Tv 


From eq. (3) and figure 5 the strain distribution is directly related to the 
local temperature T rather than the integral of T. Thus, to obtain 
accurate results for strains and stresses one would need as fine a structural 
F.E. grid as the thermal F.E. grid. Analytical results for stress from eq. 4 
give o = -179.8 MPa (-26.07 ksi); conventional F.E. results using SPAR for 2, 
4, 8, and 16 elements per length converge from -141.18 MPa (-20.475 ksi) for 
the 2-element grid to -184.38 MPa (-26.741 ksi) for the 16-element grid. A 
comparison of analytical, conventional F.E., and Ojalvo's method for stresses 
is shown in figure 5 at the right. Ojalvo's method was employed using strain 
results from two different F.E. model, the 4-element and the 8-element 
models. From the figure note that conventional F.E. results (dashed line) are 
comparable to results using Ojalvo's method (shown by symbols). In some 
instances Ojalvo's method is more accurate than conventional results and in 
some instances (where the coarse strains are not close to exact values) 
Ojalvo's method is less accurate. For other boundary conditions similar 
results were obtained; the strains varied as the local temperature and no 
benefit was realized using Ojalvo's method. 
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COMPARISON OF STRESS CALCULATION METHODS 
FOR CLAMPED ROD PROBLEM 



SPAR F.E. 



T=T Sin T = 422K(300°F) 

0 to, 

E = 72.4 GPa (10.5 x 10° psi) 
a = 23.4 X 10^ per K (13 x lo^ per °F) 

EXACT 

SPAR F.E. 

-O- • OJALVO'S IVtETHOD (4 ELEMENTS! 
-o- OJALVO'S METHOD (8 ELEMENTS) 



Figure 5 
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THIN FREE MEMBRANE SUBJECT TO A PARABOLIC 
TEMPERATURE DISTRIBUTION 


The problem shown in figure 6 is that of a thin rectangular titanium plate of 
uniform thickness, free to expand, in which the temperaure is an even function 
of y: 


T = To(l - y2/c2) 


From reference 6 


~ W 0=0 

From the thermoelastic constitutive relations and the above equations 


+ VC - (1 + v)«T] = ^ J^^aTEdy - aTE 
1-v^ ^ 

E 

o., = 


y O + ve - (1 + v)aT] = 0 

y i-v2 y y- 


solving for cx and ey gives 


= ^“^0 = Constant 


n 

£y = (l+v)aT - -jVaTo 


From the above equations the strain in the x-direction is a constant, but the 
strain in the y-direction is directly related to the local temperature T. 
Hence for accurate results for accurate values of ex, ey, and T 
are needed; since the cy varies as the local temperature in the 
y-direction as fine a discretization as that used in the thermal analysis is 
needed. If accurate or converged values of ey are not used, Ojalvo's 
method will produce less accurate results for stresses as seen in the previous 
problem. Hence if the linear strain distribution results of the 8x1 grid 
(triangular symbol) are used with the fine or exact temperature distribution 
(Ojalvo's method) the results for stresses are not as accurate as the 
conventional F.E. method in regions where strains are not converged as shown 
in the right side of the figure. Also, conventional F.E. results for the 8 x 
2 grid are close to the exact solution; therefore if Ojalvo's method is used 
with this grid, it will give little improvement. 
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THIN FREE MEMBRANE SUBJECT TO A PARABOLIC 

TEMPERATURE DISTRIBUTION (REF. 6) 

2, 2. T = 644 K(700°F) 

(l-y7c ) 0 



E = 110GPa (16xl0^psi) 

a = 10.8 
V » 0.29 


a = 10.8 X lo'^ per K (4.8 x 10 ^ per ®F) 


1 ' 

EXACT 

F.E.GRID 
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Figure 6 
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CANTILEVERED TITANIUM BOX BEAM WITH 
TEMPERATURE VARYING THROUGH THE DEPTH 


A titanium box-beam was idealized as a built-up structure consisting of 
rod and membrane elements. The fine-grid F.E. model consists of 1555 nodes, 
1512 membrane elements and 216 rod elements representing one-half of the 
entire box beam. The load applied to the structure is the temperature which 
varies in the y-direction only as shown in figure 7. This temperature 
distribution (solid curve) is a least-squares fit of a fifth-degree polynomial 
curve to the temperature distribution of reference 1. The dashed lines 
represent integrated average values of temperature which were used in an 
approximate beam theory solution to compare with F.E. results from several 
different models. 



Figure 7 




COMPARISON OF OJALVO'S METHOD WITH CONVENTIONAL F.E 
FOR A CANTILEVERED BOX BEAM 


The box beam is cantilevered as shown in figure 8. Finite-element results for 
strain in the membranes at Z = 0 at various stations along the X-axis using 
the fine grid (48 x 12 grid) are shown as the dashed curves. Beam theory 
results are shown as a solid line, and results of the coarser grids are shown 
by the symbols. Far from the ends of the beam, plane sections remain plane; 
the strains are very nearly linear and results are close to predicted beam 
theory results. At the free end (x/t = 1.0), however, the strain distribution 
is similar to the temperature distribution; this is necessary to insure the 
stress-free end condition. 

Nondimensional ized stress results at the center of the beam (x/t = 0.5) are 
shown in figure 8 at the right. Notice that results from the coarse 
structural grid (diamond symbol and dashed line) miss critical (peak) stresses 
located within the element predicted by beam theory (solid curve) and the 
fine-grid F.E. results (circular symbol). The use of Ojalvo's method, shown 
by the single-dashed curve, gives a good representation of the actual stress 
distribution consistent with results of reference 1. 




BEAM THEORY 

O SPAR F.E. (48 X 12 GRID) x/t = 0.5 
□ SPARE.E.(8x4GRID)x/l = 0.5 
0--SPAR F.EJ4x3GRID)x/t = 0.5 



Figure 8 
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COMPARISON OF OJALVO'S METHOD WITH CONVENTIONAL F.E 
FOR A CANTILEVERED BOX-BEAM (FREE END) 


A comparison of the exact solution (ox=0) with conventional F.E. and 
Ojalvo's method for the free edge of the beam (x/t=1.0) is shown in figure 9. 
Since the strain distribution at x/£=1.0 is similar to the temperature 
distribution, it is not surprising that the coarse-grid F.E. results 
for e are inaccurate. For this particular case the results from Ojalvo's 
method are worse than the results using the Conventional F.E. with a coarse 
grid. 



Figure 9 
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TEMPERATURE DISTRIBUTION IN SKIN AND STIFFENER OF A 
HEATED STIFFENED PANEL 


Another application for which Ojalvo's method might be useful, other than for 
wing and fin-type structures, would be in the stiffeners of a heated stiffened 
panel (ref. 7). As indicated in reference 8, stresses in stiffeners can be 
critical during transient heating. Shown in figure 10 are temperature 
distributions in the skin and stiffener of a stiffened panel subjected to a 
hypersonic heating simulation (ref. 7). The specimen was cooled prior to 
heating to simulate a cold soak condition which explains the sub-ambient 
temperatures over a portion of the stinger. These temperatures were used as 
input for a F.E. structural analysis. 


STIFFENER TE/VIPERATURE, °F 
600 400 200 0 -200 



STIFFENER TEMPERATURE, K 


SKIN 


600 - 
500 - 


TEMPERATURE, 

K 400 

300 


600 

400 

TEMPERATURE, 

200 


-1 0 


Figure 10 
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COMPARISON OF OJALVO'S METHOD WITH CONVENTIONAL F.E. 
FOR A STIFFENED PANEL 


A comparison of strains of the fine- and coarse-grid models indicates that the 
coarse grid is sufficient in obtaining a good approximation to the linear 
strain distribution. The coarse-grid results for stresses (solid symbol and 
dashed curve) badly miss the peak stresses which occur in the center of the 
element as shown in figure 12 by the solid line. However, when the coarse 
grid strains are coupled with the detailed temperatures (Ojalvo's method) a 
good representation of the actual stress distribution (the square symbols) is 
obtained. This indicates that stiffeners in a heated panel are appropriate 
applications of Ojalvo's method in thermal /structural analysis and design. 


RESULTS ARE THROUGH THE STIFFENER DEPTH (Z) 



Figure 12 
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CONCLUDING REMARKS 


The usefulness of a thermal/structural analysis technique for improving 
thermal stress calculations termed Ojalvo's method was investigated by 
numerical examples of several classes of thermostructural problems (fig. 13). 
The problems investigated include a rod, a thin membrane, a box beam, and a 
stiffened panel. The basis of Ojalvo's method is an observation that "strains 
in heated structures idealized by conventional components are generally less 
sensitive to spatially distributed temperature variations than are their 
corresponding stresses." Results of most bending-type problems indicate that 
Ojalvo's method is useful since the strains are related to the integrals of 
temperature and hence are less sensitive to local temperature variations. 

This means that for those problems where Ojalvo's method is appropriate the 
structural F.E. idealization may be coarser than the thermal idealization and 
also that the same coarse structural representation can be used for many 
different time slices in a quasi -static thermostructural analysis. For plane 
stress or membrane type problems the strain distributions are similar to the 
temperature distribution and a finite-element grid fine enough to calculate 
accurate temperatures would be necessary to calculate accurate stresses. This 
negates the usefulness of Ojalvo's method for this class of problems. Several 
useful areas for application of Ojalvo's method include built-up structures 
which can be idealized as bending elements (beams or plates) and stiffeners in 
a stiffened panel. 


• EVALUATE OJALVO'S METHOD FOR THERMAL STRESS ANALYSIS 

• KEY FEATURE - STRAINS ARE LESS SENSITIVE TO TEMPERATURE VARIATIONS 

THAN STRESSES 

• PROBLEMS ANALYZED 

• ROD 

• MEMBRANE 

• BOX BEAM 

• STIFFENED PANEL 

• FOR MEMBRANE AND PLANE STRESS PROBLEMS STRAINS ARE RELATED TO 
LOCAL TEMPERATURES AND OJALVO'S METHOD IS GENERALLY NOT USEFUL 

• FOR BENDING PROBLEMS STRAINS ARE RELATED TO THE INTEGRALS OF 
TEMPERATURE AND OJALVO'S METHOD IS BENEFICIAL 


Figure 13 
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